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Abstract. 

The announcement by BICEP2 of the detection of B-mode polarization consistent with pri¬ 
mordial gravitational waves with a tensor-to-scalar ratio, r = 0 . 2 ]’Iqq 5 , challenged predictions 
from most inflationary models of a lower value for r. More recent results by Planck on po¬ 
larized dust emission show that the observed tensor modes signal is compatible with pure 
foreground emission. A more significant constraint on r was then obtained by a joint analy¬ 
sis of Planck, BICEP2 and Keck Array data showing an upper limit to the tensor to scalar 
ratio r < 0.12, excluding the case r = 0 with low statistical signihcance. Forthcoming mea¬ 
surements by BICEPS, the Keck Array, and other CMB polarization experiments, open the 
possibility for making the fundamental measurement of r. Here we discuss how r sets the 
scale for models where the dark matter is created at the inflationary epoch, the generically 
called super-heavy dark matter models. We also consider the constraints on such scenarios 
given by recent data from ultrahigh energy cosmic ray observatories which set the limit on 
super-heavy dark matter particles lifetime. We discuss how super-heavy dark matter can be 
discovered by a precise measurement of r combined with future observations of ultra high 
energy cosmic rays. 
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1 Introduction 

Most of the mass in the Universe is composed of Dark Matter (DM), matter observed only 
through its gravitational interaction and not through radiation. Both astrophysics and cos¬ 
mology provide ample evidence for the existence of DM and its prevalence over visible matter 
(see e.g., [1, 2]). From the 1930’s, astrophysicists have observed dynamics in large systems 
that can only be explained by dark matter such as the dynamics of galaxy clusters [3], 
galaxy rotation curves, and more recent multi-wavelength and gravitational leasing studies 
of interacting clusters, such as the “bullet cluster” [4]. On the cosmological side, global fits 
to cosmological parameters enable a precise determination of the amount of DM present in 
the Universe. A recent result comes from Planck observations of the cosmic microwave back¬ 
ground (CMB) which find the mean DM density in the Universe in units of the critical density 
to be ^Duh'^ = 0.1198 ± 0.0015 [5], with h = 0.678 ± 0.009 kms-iMpc"\ so ^dm = 0.261. 
These indications, not only provide evidence for the presence of some unknown form of mat¬ 
ter, but also indicate that this matter should be in the form of non-relativistic (cold) particles 
[5]. 

The leading paradigm to explain DM observations is based on the Weakly Interactive 
Massive Particle (WIMP) hypotheses, which consists of two basic assumptions: (i) WIMPs 
are stable particles of mass that interact weakly with the Standard Model (SM) particles; 
(ii) WIMPs are thermal relics, i.e. they were in Local Thermal Equilibrium (LTE) in the early 
Universe. Imposing that the WIMP density today is at the observed level, using a simple 
unitarity argument for the WIMP annihilation cross section aann oc 1/M^, one obtains an 
estimate of the WIMP mass in the range of 10^ to 10^ GeV. This result, also called the 
WIMP “miracle”, links the DM problem to the new physics scale expected in the context of 
the “naturalness” argument for electroweak physics. Triggering in this way the strong hope 
that the search for WIMP DM may be connected to the discovery of new physics at the TeV 
scale. 

Searches for WIMP particles are ongoing through three different routes: direct detec¬ 
tion, indirect detection, and accelerator searches (see e.g., [2]). None of these efforts have 
discovered a clear WIMP candidate so far. In addition, no evidence for new physics has 
been observed at the Large Hadron Collider (LHC). Although not yet conclusive, the lack 
of evidence for WIMPs may imply a different solution for the DM problem outside of the 
WIMP paradigm. 

In this paper we will reconsider the scenario based on particle production due to time 
varying gravitational fields: the so-called Super Heavy Dark Matter (SHDM) scenario. This 
alternative approach is based on the possibility of particle production due to the non-adiabatic 
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expansion of the background space-time acting on the vacuum quantum fluctuations. In 
quantum theories the possibility of producing particle pairs by the effect of a strong (classical) 
external field is well known: for instance, consider the case of pair creation by strong 
electromagnetic fields. The idea to apply such a mechanism in cosmology through external 
strong gravitational fields dates back to E. Schrodinger in 1939 [6]. 

The construction of a coherent theory of particle creation by the expansion of the 
Universe was developed in the last 40 years and started with the pioneering works of [7, 8, 
9, 10, 11]. More recently, in the framework of inflationary cosmologies, it was shown that 
particle creation is a common phenomenon, not tied to any specific cosmological scenario, 
that can play a crucial role in the solution to the DM problem as SHDM (labeled by X) 
can have Dx(to) ^ 1 [12, 13, 14, 15, 16, 17, 18]. This conclusion can be drawn under three 
general hypotheses: (i) SHDM in the early Universe never reaches LTE; (ii) SHDM particles 
have mass of the order of the inflaton mass, and (hi) SHDM particles are long-living 
particles with a lifetime exceeding the age of the Universe, tx ^ to- 

Precision measurements of CMB polarization and observations of Ultra High Energy 
Cosmic Rays (UHECR) up to energies > 10^*^ eV enable a direct experimental test of the 
three SHDM hypothesis listed above. The aim of the present paper is to discuss this issue in 
light of the latest observations. 

The observation of CMB fluctuations can be directly linked to the primordial density 
perturbations in the early Universe, which provide the initial condition for structure forma¬ 
tion (see e.g. [19, 20]). Density perturbations are of two types: curvature and isocurvature, 
the first one is connected with the total energy density in the early Universe while the second 
one to the actual composition of the energy density itself. The fact that SHDM particles 
with mass Mx ~ never attain LTE implies a non negligible contribution of the parti¬ 
cles creation mechanism to isocurvature perturbations, that become source of gravitational 
potential contributing to the primordial gravitational wave background. 

Original reports from BICEP2 signalled the possibility of a substantial contribution of 
tensor modes to the CMB fluctuations [21] providing additional motivation for the SHDM 
production mechanism. As we discuss in section 2, the level of tensor modes enables the 
determination of the energy scale and mass of the inflaton held resulting in a present time 
SHDM density at the required cosmological level. The result presented by BICEP2 appears 
to be seriously polluted by the polarised light emitted by dust at galactic scales [22]. The 
recent joint analysis by BICEP2, Keck Array and Planck [23], validated on simulations of a 
dust-only modelling and performed by a simple subtraction of scaled spectra, have set the 
scale of r at r < 0.12 at 95% conhdence level. The hnal result of this analysis is expressed 
in terms of a likelihood curve for r that peaks at r = 0.05 but disfavours zero with a scarce 
statistical signihcance. As shown in [23], this result can occur by chance 8% of the time, as 
conhrmed by dust-only simulations, and cannot be considered as a conclusive detection of 
primordial B-modes. In the near future several different detectors, both ground-based and 
sub-orbital, will take data at a variety of frequency enabling a more precise determination of 
the CMB primordial tensor modes. 

The third and last hypothesis of the SHDM paradigm is a long particle life-time, much 
longer than the age of the Universe. This is a general requirement of any DM model based 
on the existence of new particles, common to the WIMP hypothesis as well. As in the case 
of WIMPs, discrete gauge symmetries protecting particles from fast decays work equally well 
for SHDM particles. A general particle physics example of SHDM with life-time exceeding 
the age of the Universe can be found in [24]. 
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The best way to test the existence of SHDM is through the indirect detection of its 
decay and/or annihilation products since direct detection of SHDM is unattainable. For 
instance, the annihilation of SHDM can occur in the case they are ”majorana-hke” particles, 
i.e. particles and anti-particles coincide. Since the annihilation cross section of a (point) 
particle is bounded by unitarity, aann oc 1/Mi 1/M^, the annihilation process results in 
an unobservably small annihilation rate [25]. Even if alternative theoretical models can be 
constructed, with a sizable annihilation rate [26], we will not consider this case here and 
focus on the more general case of signatures from SHDM decays. 

If SHDM particles decay, we can determine the composition and spectra of the SM 
particles produced by SHDM decays under general assumptions on the underlying theory 
[25, 27]. Given the mass Mx ~ and the lifetime tx to, UHECR experiments with 
high statistics should be able to detect the decay products at extreme energies > 10^*^ eV. In 
section 3 we discuss the possibility of indirect detection of SHDM particles focusing on spectra 
and chemical composition taking into account the latest UHECR observations performed by 
Auger [28, 29, 30] and Telescope Array (TA) [31]. 

When CMB polarization experiments determine the inflation scale, UHECR detectors 
will study the parameter space for SHDM lifetimes. In section 4 we will draw our conclu¬ 
sions discussing the future possibilities in constraining SHDM models through CMB and 
UHECR observations, highlighting the discovery capabilities of the next generation UHECR 
experiments. 

2 Tensor modes and Super Heavy Dark Matter 

In the scenario based on gravitational production, SHDM never attains LTE and we can 
connect its density today with the density when it was created. Following [12] we can write 

Pxjtp) _ PxjtRH) ^2 1 ) 

PR{to) PR{tRH) \ To J 

where pn is the energy density of radiation, px is the SHDM energy density, and Trh and 
To are the temperature of the Universe at reheating time tun and today to, respectively. 
Assuming that SHDM particles are produced just after the de Sitter phase at time tf., when 
the coherent oscillation phase begins, then the inflaton energy density and the SHDM energy 
density will be redshifted at approximately the same rate. This holds until the reheating 
phase hnishes and the radiation dominated phase begins. Therefore, assuming that before 
reheating most of the energy density of the Universe is converted into radiation [12] (p/j ~ 
Pc = SH'^Mpi/Sir), we can rewrite 

PxjtRH) Stt Pxjte) /2 21 

PRjtRn) 3 

where Mpi = 10^® GeV (the Planck mass), and Hjtf.) is the Hubble parameter at tg- 

The energy density of particles produced by the interaction of a (time varying) classical 
gravitational field with a quantum vacuum depends on the details of the underlying quantum 
theory (e.g., the choice of vacuum) and cosmology (e.g., the structure of the time dependence 
of the scale factor). These details were analyzed by several authors in the framework of 
different models, however, as was shown in [12, 13, 32], the SHDM density at the end of 
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inflation can be generally approximated by [32, 33] 

«(*.) ~ , (2.3) 

quite independent of the specific underlying theory. Eq. (2.3) is inaccurate in the low mass 
limit Mx 'C M^j,, while it works reasonably well near the cut-off region Mx ^ that is 
the regime that we are interested in (see later). For a more accurate determination of the 
SHDM density, obtained by numerically solving mode equations, taking into account also the 
isocurvature-curvature cross correlation, see [34, 35]. 

Using Eq. (2.1), (2.2) (2.3), and taking the Hubble parameter at equal to the inflaton 
mass H{tg) ~ one gets a simple estimate of the ratio of the SHDM density to the critical 
density today Qx = Px{to)/ Pc{to)- 

being Tq = 2.3 x 10“^^ GeV the CMB temperature today and Qr = 4 x 10“® the radiation 
density today. 

The gravitational production of SHDM is intimately connected to the process that, dur¬ 
ing inflation, generated primordial large-scale density fluctuations. As realized in [33], the 
observations of the CMB fluctuations provide important insights into SHDM properties hav¬ 
ing the potential to falsify the model. The production of SHDM during inflation gives rise to 
isocurvature perturbations that become sources of gravitational potential energy contributing 
to the tensor power spectrum of the CMB. 

The observation of a sizeable amount of tensor modes in the CMB fluctuation pattern 
would imply an important confirmation of the inflationary paradigm in agreement with the 
hypothesis that adiabatic perturbations originate within the single field, slow-roll framework 
of inflation. 

To asses the impact of primordial B-modes on SHDM, we express the inflaton mass 
in terms of the tensor to scalar ratio r. At leading order in the slow roll approximation we 
can write the inflaton potential height 4* in terms of r and the amplitude Ag of the scalar 
perturbations on super-Hubble scales. Using the combined analysis of Planck and WMAP 
data one has [36] : 

U* ~ ~ (2 X 10i6QgV)4^ = M^uT (2.5) 

being Mgut = 2 x 10^® GeV the scale of Grand Unification Theory (GUT) and ro = 0.12 
is the experimental upper bound on the tensor to scalar ratio [36]. It is interesting to note 
here that in the case of a substantial contribution of B-modes, as in the case of the BICEP2 
claim [21], one finds an inflaton potential height intriguingly near the typical GUT scale. Let 
us now assume an inflaton potential of the type: 

( 2 . 6 ) 

where /3 = 2/3,1,4/3, 2. From the latest Planck observations [36] follows that the quadratic 
potential is moderately disfavoured while higher values of (5 are strongly disfavoured. Values 


- 4 - 









Figure 1. [Left panel] Inflaton mass as function of the ratio r of tensor to scalar modes for different 
choices of the inflaton potential as labeled. [Right panel] Ratio of the SHDM mass and inflaton mass 
as function of r, obtained as solution of the equation fix = ^dm using different choices of the inflaton 
potential as labelled. 


of /? = 2/3,1,4/3 are motivated by axion monodromy, which combines chaotic inflation and 
natural inflation [37, 38, 39, 40]. 

From Eqs. (2.5) and (2.6), using the definition of the slow-roll parameter e((/>) 


e(</) 


Mh 

IGvr V{cj)) 



(2.7) 


we can determine the inflaton mass as a function of the tensor to scalar ratio r: 


= Mgut 


/ Mgut Y ( (— 

Mpi ) V ^ ) V^o 


\ 1+/5/2' 


4-/J 


( 2 . 8 ) 


In figure 1, left panel, we plot the value of the inflaton mass as function of r, choosing 
different values of the inflaton potential power law index /? = 2/3,1,4/3, 2 [36, 37, 38, 39, 40]. 

Using Eq. (2.8) we can rewrite the density of SHDM today in terms of the inflaton 
potential power law index /3, the ratio of tensor to scalar modes r and the ratio Mx/M^-. 

(ISK'j X 


8 _ 2/3 2+/3 

2 / Mgut \ f \ / r\ 4-/3 

[ MpiJ [ /3 J [roj 


mJ 


(2.9) 


Taking the reheating temperature TpH — 10® GeV and assuming that the SHDM density 
today coincides with the observed DM density Dx = ^dm = 0.261, using Eq. (2.9), we can 
determine the ratio Mx/M^p as function of r for different choices of the inflaton potential. 
In the right panel of figure 1 we plot the resulting behaviours. As already discussed in [33], 
although only in the case of an inflaton potential with /3 = 2, the typical values of the SHDM 
mass are such that Mx/M^ > 0(5) with a range in mass that spans from 10^^ GeV up to 
10^^ GeV, depending on the value of r and the choice of the inflaton potential power law 
index (3. 
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The result presented in the right panel of figure 1 also depends on the assumption about 
the reheating temperature. For instance, in the case /3 = 2 the equation Qx = ^dm = 0.261 
has no solutions for r < 3 x 10“^, in this case a larger value^ of Tjih, at the level of Tjih — 10^*^ 
GeV, is needed in order to obtain the correct SHDM density today for r < 3 x 10“^. 

To assess the dependence on the reheating temperature of the result presented in the 
right panel of figure 1, we have repeated that computation using different values of Trh- In 
the case /3 = 2, as discussed in [33], there is a lower bound on Trh, fixing Trh at its minimum 
the equation fix = n_DM shows no solution only for a tensor to scalar ratio r < 2 x 10“^, 
while at larger r the ratio Mx/M^ differs by no more than 50% from the result of figure 1. 
In the case /3 = 2/3,1,4/3 there is no lower bound on the T^h value. In these cases, to check 
the dependence of our results on T^h we have repeated the computations of figure 1 (right 
panel) taking different values of Trh down to 10® GeV always obtaining solutions Mx/M^ 
that differ by no more than 50% from the values of figure 1 (where Tjih = 10® GeV). A 
change in the Mx parameter by a factor of two, changing the Trh value, does not affect the 
result presented here, therefore in the following we will consider only the case Tjih = 10® 
GeV as a reference value. 

Another important point concerns the limits on the the SHDM models that come from 
the experimental limits on isocurvature perturbations. Almost independently of Tjih, we 
found that the ratio Mx/M^ needed to obtain the observed DM density is in the range 2-^15. 
With this parameters choice, as shown by [33] and more recently by [34], the isocurvature 
perturbations expected are well within the experimental limits. 

We can conclude this section by stating that super-heavy particles produced by time 
varying gravitational fields in the early Universe provide a viable explanation of the DM 
problem even with a ratio of tensor to scalar modes at the level of 10“®. This is a non 
trivial result that shows how the observation of a non zero fraction of tensor modes in the 
CMB fluctuations pattern, already at the level of confirms that the gravitational 

production of SHDM particles is a viable mechanism to explain the DM problem, assuring a 
density of SHDM today at the observed level. In the next section we will consider the case 
of SHDM discussing its possible indirect detection through UHECR observations. 

3 Ultra High Energy particles by Super Heavy Dark Matter decay 

In order to constitute a substantial fraction of the observed DM, SHDM should be composed 
by metastable particles with lifetime longer or much longer than the age of the Universe. In 
the previous section we have discussed the gravitational production of SHDM showing how 
the mass of these particles is fixed by cosmology. We are left with a single free parameter 
in the model: the particle lifetime rx, that can be constrained through observations of the 
decay products. 

Following [25, 27], on very general grounds, we can assume that the SHDM decay 
produces a pair of quark anti-quark that gives rise to a parton cascade, producing SM particles 
through the subsequent hadronization process. The basic signatures of the SHDM decay 
process are three: (i) SHDM particles (as any other DM particle) cluster gravitationally 
and accumulate in the halo of our Galaxy with an average density of ~ 0.3 GeV/cm®; 
(ii) in the hadronic cascades the most abundant particle produced are pions, therefore UHE 

^In inflationary scenarios embedded in gravity mediated SUSY breaking the gravitino over-production 
bound restricts the reheating temperature to Trh ;$ x 10® -h 10^® GeV that in the case (3 = 2 corresponds to 
a lower bound on the tensor to scalar ratio around 10“® (see [33] and references therein). 
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neutrinos and photons are the most abundant particles expected on Earth; (iii) the non¬ 
central position of the sun in the Galactic halo results in an anisotropic flux of the SHDM 
decay products [41, 42], 

The quantitative predictions for energy spectra and chemical composition of UHECR 
coming from SHDM decay require an extrapolation of QCD calculations from the TeV scale 
up to the scale Mx- There are several computational schemes proposed in the literature based 
on analytic approximations [43] or Monte Carlo simulations [27, 44, 45, 46, 47]. These schemes 
predict accurately the secondary spectra of SM particles produced in the SHDM decay and 
agree well each other. Their most important outcome from the observational point of view 
is a flat spectrum, that at the relevant energies can be approximated as dN/dE oc 
independently of the particle type, with a photon/nucleon ratio of about 'j/N ~ 2 3 and a 

neutrino/nucleon ratio u/N ~ 3 -i- 4, quite independent of the energy [27]. 

The UHECR emissivity produced by the decay of SHDM in the halo of our Galaxy can 
be written as 


Ip,'y,u{E) — 


1 dNp^^^^{E) 

Mxtx dE 


Px{R) 


(3.1) 


where Mx and tx are the mass and lifetime of the SHDM particle, dN/dE is the energy 
spectrum (fragmentation function) of the SM particles produced by the decay and pxiR) 
is the density of X-particles in the Galaxy as function of the distance R from the Galactic 
Center (GC). 

In the present paper, assuming that SHDM contributes to a substantial fraction of the 
galactic DM, we will use the DM density profile as computed in several numerical simulations 
by Navarro-Frenk-White (NEW) [48] and Moore [49]: 


(R/R,)“(l + R/R,)3-“ 

with OL = 1,3/2 in the case of NEW and Moore respectively and Rg = 45 kpc being the 
fiducial radius of the DM distribution as obtained in [50]. The DM density is normalized to 
the expected local (solar) density, namely: px{Rq) = P©^ = 0.3 GeV/cm^. 

The expected flux on Earth is given by the integral over the line of sight of the emissivity 
(3.1), indicating with 6 the angle between the line of sight and the axis Earth-GC we can 
write the flux as 


with 


1 

,R{E,e) = —j^ dsIi{E,s) 
Smax{0) = Rq cos{6) + R% + R^siii^O, 


(3.3) 


being i = p,'y,iy the type of particle produced in the SHDM decay, s the coordinate along the 
line of sight, Rq = 8.5 kpc the distance between Earth and GC and Rh = 100 kpc the radius 
of the galactic halo. Changing integration variable in equation (3.3) from the line of sight 
coordinate s to the galactic radius R, one obtains a more convenient formula for numerical 
computations [25]: 


ME, 6) 


1 dNijE) 
AtttxMx dE 


jRQsinie) J- R%siTi?{e) 


(3.4) 
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Together with UHECR emitted by the decay of SHDM in our local halo (3.3) one 
should take into account also the contribution of the overall Universe to the emitted flux 
of particles. In the case of protons and photons the energy losses due to the interactions 
with astrophysical backgrounds (mainly CMB) reduce the expected cosmological flux to a 
negligible fraction with respect to the flux produced by our local halo [25]. In the case of 
neutrinos the contribution to the expected flux from the whole Universe is not negligible 
being at a level of about 10% of the local flux. Therefore one should add to the neutrino flux 
given by (3.3) also the flux given by: 




C ^XPc 

Att Mxtx 



dN^ 

dE 


[(1 + 


(3.5) 


where dt/dz = —HqsJ (1 + z)^Qrn + Pc = 5.5 x 10“® GeV/cm^ is the critical density of 
the Universe and Su{Ey,z) is the neutrino opacity of the Universe calculated in [51, 52]. 

The contribution of SHDM decay to the UHECR fluxes starts to be relevant at the 
highest energies (E > 5 x 10^® eV). Using the model proposed in [53] for UHECR by astro- 
physical sources (see also [54, 55, 56, 57, 58, 59] for different computation schemes and models 
of UHECR composition), in hgure 2 we have plotted the total UHECR flux, highlighting the 
contribution of the different components: proton (dotted green), He (dotted magenta), CNO 
(dotted blue), MgAlSi (dotted cyan), Fe (dotted black) from astrophysical sources [53] and 
proton (green solid), photon (blue solid) and neutrino (red solid) from the decay of SHDM. 
The latter fluxes where normalized integrating over the whole sky (0 ^ 0 < tt), with a Moore 
density prohle for SHDM. The four panels of hgure 2 correspond to the four different as¬ 
sumptions on the inhaton potential discussed in the previous section [3 = 2/3,1,4/3, 2 (as 
labelled in the hgure) and hxing a reference value of the tensor to scalar ratio r = 0.05, 
that corresponds to the peak in the tensor to scalar ratio likelihood curve of the combined 
analysis of Planck, BICEP2 and Keck Array [23]. As labelled in the hgure, this choice of 
the inhation parameters corresponds to a SHDM mass: Mx = 4.5 x 10^^ GeV in the case 

of ^ = 2, Mx = 4.3 X 10^5 GeV with /3 = 4/3, Mx = 1.7 x lO^^ Q^y = 1 and 

Mx = 4.7 X 10^® GeV with /3 = 2/3. The minimum lifetime of SHDM was hxed by imposing 
the Auger experimental limits on the 7 ray hux at E > 10^® eV, observed at the level of 
2% [60, 61], and on the number of events at E > 6 x 10^® eV [28]. These minimum lifetime 
values are: tx = 2.2 x 10^^ yr in the case of /3 = 2, rx = 1.2 x 10 ^^ yr in the case of /I = 4/3, 

/3 = 8.3 X 10^^ yr in the case of /3 = 1 and tx = 6.3 x 10^^ yr in the case of /3 = 2/3. 

Let us now take into account the composition of UHECR. The chemical composition is 
inferred by the mean value of the depth of shower maximum (Xmax) and its dispersion (RMS) 
(^{Xmax)- As was hrst discussed in [62] (see also [63]), the combined analysis of {Xmax) and 
<j{Xmax) is more sensitive to the chemical composition and provides less model dependent 
results. The main uncertainties in such a procedure are provided by the dependence of {Xmax) 
and its fluctuations on the interaction models used for the shower development. Most of 
such models ht low energy accelerator data while providing somewhat different results when 
extrapolated to the energies of relevance for UHECRs (for a review see [64] and references 
therein). 

In our calculations we follow the procedure of [65] where four models of HE interaction 
were included to describe the atmospheric shower development, namely EPOS 1.99 [66], 
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Figure 2. UHECR flux: proton (dotted green), He (dotted magenta), CNO (dotted blue), MgAlSi 
(dotted cyan), Fe (dotted black) from astrophysical sources [53] and proton (green solid), photon (blue 
solid) and neutrino (red solid) from the decay of SHDM with a Moore density profile. The total flux 
is represented by the black thick solid line. Experimental data are the latest observations of Auger 
[28] and TA [31]. All plots are obtained assuming r = 0.05, taking the four different choices of the 
inflaton potential: (3 = 2 upper left panel, P = 4/3 upper right panel, /3 = 1 lower left panel and 
/? = 2/3 lower right panel. The corresponding values of the SHDM parameters {Mx,tx) are labelled 
in the different panels. 


Sibyll 2.1 [67], QGSJet 01 [68] and QGSJet 02 [69], in order to derive for each given nuclear 
primary a simple prescription for {X^ax) and a{Xmax)- To account also for the contribution 
of UHE photons to X^ax and cr{Xmax) we have used GONEX simulations of y-induced 
showers, taking into account the Landau Pomeranchuk Migdal (LPM) and the Geomagnetic 
field effects on the showers development [70]. 

The observations of Auger point toward a mixed composition of UHEGR with a preva¬ 
lent light composition at low energies E < 5x 10^® eV and a progressively heavier composition 
at the highest energies [28, 29, 30]. The highest energy bin with an observed chemical compo¬ 
sition is placed at E ~ 5 x 10^® eV, at highest energies there are no available measurements 
of composition [28, 29, 30]. As discussed above, the decay of SHDM gives a substantial 
contribution to the composition of UHECR only at energies larger than 5 x 10^® eV, where 
the fraction of photons and neutrinos rises over a few %, therefore in an energy band not 
scrutinized by Auger and TA because of the extremely low statistics of the collected events. 

In hgure 3 we plot (Xmax) and cr{Xmax) corresponding to the fluxes of the upper left 
panel of figure 2, obtained with a tensor to scalar ratio r = 0.05 and an inflaton potential 
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Figure 3. {Xmax) (left panel) and cr{Xmax) (right panel) as computed in this paper compared 
with the latest Auger observations [28, 29, 30]. Shown by the red band is the uncertainty in the 
computations due to the choice of the hadron interaction model (see text). 


power law index /3 = 2, together with the Auger data of 2013 (black points) [28, 29, 30]. Apart 
from the last energy bin in the RMS data, we can conclude that a mixture of astrophysical 
sources and a component produced by the decay of SHDM is compatible with the chemical 
composition observed by Auger. The chemical composition expected with different choices 
of /3 remains practically unchanged, therefore in figure 3 we have plotted only the case /3 = 2. 

Future observatories of UHECRs and neutrinos should be able to discover SHDM or 
constrain their lifetimes. For UHECRs, the Extreme Universe Space Observatory to be 
deployed in the Japanese Experiment Module of the International Space Station (JEM- 
EUSO) should achieve about an order of magnitude higher exposure at the 10^® eV energy 
region as compared to current ground arrays [71]. A possible future network of such space 
observatories could continue to improve sensitivity up to the observational limit of the Earth’s 
atmosphere, which is another four orders of magnitudes in exposure. For UHE neutrinos, the 
upcoming Askaryan Radio Array (ARA) [72, 73] will improve by about an order of magnitude 
the sensitivity to 10^^ eV when compared to the strongest current limits of ANITA H [74, 75]. 

The most striking signature of SHDM models is represented by the peculiar composition 
of UHECR at the highest energies, with a flux dominated by neutrinos and photons, which 
is an exact outcome of the sole decaying dynamics of the SHDM particles [27]. In figure 4 
(left panel) we plot the fraction of protons (only those produced by SHDM), photons and 
neutrinos over the total flux (taking again the case /3 = 2, r = 0.05 in both panels of the 
figure). 

In figure 4 (right panel) we plot the flux of UHE cosmic rays, gamma-rays, and neutri¬ 
nos from SHDM decay and the sensitivity of ARA to neutrinos (thick blue solid line) and 
JEM-EUSO (thick red solid line) to cosmic and gamma-rays. Right panel of figure 4 shows 
how neutrino observatories are not competitive with UHECR observatories in the search for 
SHDM. In the near future, only the JEM-EUSO experiment should be able to detect these 
secondaries and test the SHDM hypothesis if the predicted ratio of protons, neutrinos, and 
photons is observed (figure 4 left panel). If none of the secondaries are found, a lower limit 
on the lifetime of SHDM particles can be derived. 

In order to better quantify the JEM-EUSO capabilities in exploring SHDM models, 
in figure 5 we plot the regions on the plane {r,Tx) accessible to this experiment. While 
the SHDM mass Mx is univocally fixed by r (and the inflaton potential), to determine the 
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Figure 4. [Left panel] Fraction over the total UHECR flux of protons, photons and neutrinos by 
SHDM decay, in the case of r = 0.05 and /? = 2. [Right panel] Sensitivity to SHDM decay products: 
for UHECRs the future JEM-EUSO space mission (thick red solid line) and for UHE neutrinos the 
upcoming ARA observatory (thick blue solid line). Experimental data are those of Auger (red points) 
[28] and TA (black points) [31]. Theoretical fluxes are the same of upper left panel in figure 2 (r = 0.05, 
/3 = 2). 

detectable values of tx we assumed a JEM-EUSO exposure that enables the observation of 
UHECR fluxes from SHDM a factor of 30 lower than the value determined by the Auger 
limits on 7 rays at E > 10^® eV and on the total number of events observed by Auger at 
E > 6 X 10^® eV. The latter being the minimum allowed value of tx used in figures 2, 3 
and 4. In figure 5 we have considered the tensor to scalar ratio in the range (10“^, 10~^), as 
follows from the latest combined analysis of BICEP2, Planck and Keck Array [23]. Only in 
the case /3 = 2 we restricted our computations to the interval 3 x 10“^ < r < 10“^ because 
of our assumptions on the reheating temperature as discussed in the previous section (see 
right panel of figure 1). The four panels of figure 5 refer to the four different choices of 
the inflaton potential power law index /3 = 2/3,1,4/3,2 (as labelled in the figure). In each 
panel the lower abscissa represents the tensor to scale ratio r while the upper one shows the 
corresponding value of the SHDM mass Mx that assures Dx = D dm today. 

4 Conclusions 

The BICEP program represents great progress in the search for primordial B-mode polariza¬ 
tion. This encouraging trend revived the notion that the signal may be higher than previously 
expected. However, Plack showed that subtracting foreground contamination by dust is also 
more challenging than past estimates. The joint Planck, BICEP2, and Keck Array analysis 
[23] set an upper limit on r while mildly pointing towards a non-zero possibility. Given the 
large interest in the community and the ability of next generation experiments to surpass 
current challenges, we are optimistic that the fundamental measurement of r will be reached 
in the near future, if it is not too dehantly small. 

Relatively large values of r, which may be detected in the near future, motivate the 
idea that dark matter is mainly composed of superheavy relics from the inflationary epoch. 
Given a measurement of r, the existence of SHDM can best be tested by exploring the 
decay lifetime parameter range with future UHECR observatories. The higher the statistics 
of UHECR experiments at energies around 10^^ eV, the more likely the detection of SM 
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Figure 5. Shadowed red areas represent the regions in the plane {r,Tx) (or {Mx,tx) as labelled in 
the upper abscissa of each panel) accessible to the JEM-EUSO experiment, each panel corresponds 
to a different choice of the power law index of the inflaton potential (as labelled): upper left panel 
/3 = 2, upper right panel /3 = 4/3, lower left panel /3 = 1 and lower right panel /3 = 2/3. 


particles produced by SHDM decay. Given the currently planned UHECR and UHE neutrino 
detectors, JEM-EUSO is best placed to effectively study the allowed SHDM lifetimes with 
possible lifetime detections or constraints reaching values as high as tx — 10^^ yr. 
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